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ABSTRACT 

In many astrophysical situations, as in the coalescence of supermassive black hole pairs at gas rich galactic 
nuclei, the dynamical friction experienced by an object is a combination of its own wake as well as the wakes of 
its companions. Using a semi-analytic approach, we investigate the composite wake due to, and the resulting 
drag forces on, double perturbers that are placed at the opposite sides of the orbital center and move on a 
circular orbit in a uniform gaseous medium. The circular orbit makes the wake of each perturber asymmetric, 
creating an overdense tail at the trailing side. The tail not only drags the perturber backward but it also exerts 
a positive torque on the companion. For equal-mass perturbers, the positive torque created by the companion 
wake is, on average, a fraction ~ 40 - 50% of the negative torque created by its own wake, but this fraction 
may be even larger for perturbers moving subsonically. This suggests that the orbital decay of a perturber in a 
double system, especially in the subsonic regime, can take considerably longer than in isolation. We provide 
the fitting formulae for the forces due to the companion wake and discuss our results in light of recent numerical 
simulations for mergers of binary black holes. 
Subject headings: binaries : general — black hole physics — hydrodynamics — waves 



1. INTRODUCTION 

Understanding the nature of the dynamical friction (DF) 
force is of great importance to describe the evolution 
of gravitational syst ems. Althoug h the concept of DF 
was introduced by IChandrasekhaJ (Il943h studying col- 
lisionless backgrounds, it is also of astrophysica l in- 
terest wh en considering gaseous media (e.g, DokuchaevI 
19641: iRuderman & Spiesel 1971; RephaeH & Sabeter 1980; 



OstrikejfT999E Sanchez- Salcedo & Brandenburg 1999, 2001). 



In a seminal paper, Ostriker (1999) derived the analytic for- 
mulae for drag forces on a perturber moving straight in a 
uniform gaseous medium, which have found a vari ety of as- 
trophysical applications, from accretion disks (e.g. iNarayanI 
EOOO"; "Karas & Subr '200?; "Chang" '200?) to the in tracluster 
medium (e.g., Kim et al. 2005; Kim 2007; Conroy & Ostriker 
|2008). Recently, Kim & Kim (2007, hereafter iPagerl ex- 
tended the work of .Ostriker. (.1999) to a more real istic case 
where the perturber moves on a circular orbit (see iBaraussd 
|2()07 for the relativistic case). 

Although the works mentioned above have improved our 
understanding on DF force on a single object, there are 
many astronomical situations involving double or multiple 
bodies in which one needs an analytical estimate of the net 
drag force. While the orbital evolution of a binary sys- 
tem due to the DF in a collisionless system has been stud- 
ied in great detail (e.g., lHeggidlT975|) . it is still lacking for 
the gaseous case. The latter is crucial to descri be the for- 
matio n and hardening of close binary systems (iBate et al.l 
l2002h. me rging of double black holes at galactic centers 
(lEscala et al. 2004, 2005; Dotti et al. 2006, 2007; Mayer et al. 
l2007h . and orbital decays of kpc-sized giant clumps formed 
in primordial ''clump cluster" galaxies ( Immeli et al.. .2004 : 
iBournaud et al.ll2007l) . In particular, angular momentum loss 
to gas provides a plausible mechanism to explain the rel- 
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atively rapid coalescence of supermassive black hole pairs 
in galaxy centers (e.g., Begelman et al. 1980; Gould & RixJ 
2000; Armi tage & Natarajanl2QQ2ll2005h . lEscala et al.1 (l2004 
■2005.) found, through numerical simulations, that the black 
holes produce a composite wake of an inclined ellipsoidal 
shape that exerts a net torque on the black holes. This clearly 
indicates that the in- spiral of one black hole is affected also 
by the wake from its companion. In this Letter, we consider 
a system composed of two perturbers on coplanar circular or- 
bits in order to assess quantitatively the effect of the com- 
panion wake. Using a semi-analytic approach, we study the 
structure of the combined density wake, evaluate the resulting 
drag force on each object, and apply our results to the cases 
considered in numerical simulations of DF-induced mergers 
of black holes. 

2. FORMULATION 

We consider two point-mass perturbers moving on copla- 
nar circular orbits around the same orbital center in an invis- 
cid gaseous medium.^ In order to isolate the effects of the 
companion, we assume a perfectly uniform background with 
density po and ignore the complications from density gradi- 
ents. To simplify the presentation further, we consider that 
both perturbers orbit with a fixed radius Rp and at a constant 
velocity Vp, and are located at the opposite sides of the orbital 
center, but the extension to general case is straightforward. 
To study the response of gas to the perturbers, we employ the 
same formalism as in Paper I for a time-dependent linear per- 
turbation analysis. The reader is referred to Paper I for a more 
detailed description of the method. 

Assuming that the perturbed density a = (p-po)/po is adi- 
abatic and very small, we linearize the equations of hydrody- 
namics to obtain a three-dimensional wave equation 



V^Qf- 



1 d^a __47rG 
c^W~ 



-p^^ti^.t), 



(1) 



^ The perturbers may orbit under a common external gravitational poten- 
tial and/or under their mutual gravity. They constitute a binary if the latter 
dominates. 
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where Cg is the adiabatic speed of sound in the unperturbed 
medi um and pext de notes the mass density of the perturbers 
re.g. jOstrikei]ll999h . 

We work in cyHndrical coordinates (R^ip^z) with its origin 
lying at the orbital center and the z-axis perpendicular to the 
orbital plane. Assuming that the perturbers with mass Mp and 
fMp each (with / denoting the mass ratio) are introduced at 
(/?p, 0, 0) and (Rp^n^ 0) at ^ = 0, respectively, one can write 



p,xt(x,t) = Mpn(t)5(R-Rp)S(z) 

X {S[Rp(^-nt)]^fS[Rp(^-7r-nt)]}, 



(2) 



where H(t) is the Heaviside step function and ft = Vn/Rp is 
the angular speed of the perturbers. Since equation ([IJ) is lin- 
ear, a is given by a simple superposition of the wakes of both 
perturbers. By solving equation ([T]) based on the retarded 
Green's function technique and simplifying the resulting in- 
tegral analytically, one can show that the perturbed density is 
reduced to 



a(x,0- 



GMp 

C^Rn 



V{R,if,z,t), 



(3) 



where V = Vi+fV2 with Vi and V2 denoting the dimension- 
less wake, given by equation (8) in Paper I, of the perturber 
with mass Mp and fMp, respectively. Note that V2{R^ ^^ z, t) - 
V\ {R^ip-Ti^z^ t) for the perturbers in consideration. 

The gravitational drag force exerted on the perturber of 
mass Mp located at the position Xp can be obtained by directly 
evaluating the integral 



Fdf = GMpPo / d'x 



I' 



3 a(x,t)(x-Xp) 



^DF,l 



+ Fi 



DF,2, 



(4) 



where Fdf,i = -J^(Xi,/?R+Xi,^(^) and Fdf,2 = -J^/(X2,/?R + 
T2^ip(p), with T = 47Tpo(GMp/Vp)^. Here, R represents the 
unitary direction vector along R, and Xi^r and Xi^^ are the di- 
mensionless drag forces on a perturber by its own wake in 
the radial and azimuthal directions, respectively, while 22,/? 
and X2^cp refer to those from the wake of its companion with 
mass fMp. The former was evaluated and widely discussed 
in iPaper 1 Wq here focus on X2,/? and X2,(p defined by equa- 
tions (12) of iPaper I except T>2 replacing V for the perturbed 
density. 

As in Paper I, we calculate V(x^t) and Fdf on a three- 
dimensional Cartesian mesh centered at the center of the orbit. 
We checked that the grid spacing of ~ 7^^/640 and the box 
size of ~ (20- 100) Rp are sufficient to give converged results. 

3. RESULTS 

3.1. Density Wake 

In this section, we limit our presentation to the cases with 
equal-mass perturbers; the cases with / ^ 1 will be briefly dis- 
cussed in iH As the perturbers introduced at ^ = move along 
a circular orbit, they continuously launch sound waves that 
propagate and affect the surrounding medium that would oth- 
erwise be uniform and static. Any location inside the causal 
region is able to receive sonic perturbations from both per- 
turbers, possibly multiple times, creating a density wake that 
differs significantly depending on the Mach number Ai = 
Vp/cs of the perturbers. Figure [T] displays the distributions of 
the dimensionless wake V on the orbital plane (z = 0), when a 
steady state is reached (t > Rp/cs) for M = 0.6, 1.2, 2.0, and 
3.2. In each frame, the perturbers located at (x^y) = (±Rp^O) 
are rotating counterclockwise. 




4-4 4 

Fig. 1 . — Distributions of the dimensionless perturbed density P in a steady 
state for (a) M = 0.6, (b) 1.2, (c) 2.0, and (d) 3.2 on the orbital plane (z = 0). 
The perturbers of equal mass located at (x, _y) = (±^p, 0) are moving in the 
counterclockwise direction along a circular orbit marked by the black circle 
in each frame. Colorbar labels log P. 



For subsonic perturbers (M < 1), the perturbed density is 
smooth without involving a shock. The bending of the wakes 
caused by the circular motions leads to slight over-densities at 
the trailing sides (see Fig. [1^), producing nonvanishing drag 
forces. In the steady state, which is achieved at t ^ 00, a 
parcel of gas at any position receives one sonic perturbation 
from each perturber. For supersonic cases (M > 1), on the 
other hand, the wake of each perturber initially consists of a 
sonic sphere and a Mach cone, the interiors of which are in- 
fluenced by sonic disturbances once and twice, respectively. 
Because of the circular motion, a perturber (and the head of its 
Mach cone) is able to overtake its own sonic sphere and sub- 
sequently the other sonic sphere from the companion, both 
of which are expanding radially outward. This in turn pro- 
vides additional perturbations to the wakes and thus forms 
long high-density tails that loosely wrap the perturbers in a 
trailing spiral fashion, as Figure [l] shows. The tails in fact 
trace the regions bounded by shock discontinuities where the 
gas has received sonic signals four times (three from one per- 
turber and one from the other) and do not overlap with each 
other, provided the Mach number is less than 2.972 (see be- 
low). Note that the densest parts of a tail are located at the 
immediate trailing side of a perturber, which indicates that the 
companion wake generally tends to reduce the net DF force, 
as we shall show in ^3.21 

The wake tails thicken as the Mach number increases from 
unity. Paper I showed that a tail from a single perturber be- 
comes fat enough to make the inner edge contact with the 
outer edge at Mi= 4.603. The self-overlapping of a tail de- 
velops a new thin tail that becomes thicker with increasing A4 
and again overlaps itself at M2 = 7.790, and so on. In double- 
perturber cases, however, one tail is able to mutually over- 
lap with the other even before it undergoes self-overlapping. 
The critical Mach numbers Mn for the m utual ove rlapping 
of tails are determined by equation (B3) in iPaper J for half- 
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integer n. A few critical Mach numbers are A^i/2 = 2.972, 
M.3/2 = 6.202, and M5/2 = 9.37 L The high-density narrow 
tails shown in Figure[T]l for M = 3.2 are constructed by com- 
bining six sonic disturbances emitted by the perturbers. 



3.2. Gravitational Drag Force 

Once the perturbed density V(x^t) is constructed, it is 
straightforward to calculate the DF forces exerted on each 
perturber. While the volume of space influenced by the per- 
turbers steadily increases with time, the resulting drag forces 
quickly converge to their steady- state values typically within 
one orbital period. Figure [2l plots the various drag forces in 
a steady state on a perturber as functions of the Mach num- 
ber. In order to avoid a divergence of the force integral, only 
the region with the distance r > rmm = Rp/ 10 from the per- 
turber is taken into account in the force computation. Note 
that only Xi^^ depends on the Coulomb logarithm ln(rj^in/Rp) 
for supersonic perturbers; the other three forces (Ti^r, T2,r, 
and 22,(p) are independent of the adopted value for rmin. The 
local bumps in 22,/? and X2,(p ^^^ caused by the overlapping of 
the tails occurring at the critical Mach numbers, as discussed 
ingm 

Figure [2] shows that 22,^^ has opposite sign to Xi,^^ for all 
M . This implies that, regardless of the Mach number, one 
perturber in a double-perturber system gains angular momen- 
tum from the gravitational torque exerted by the companion 
wake, while its own wake always takes away angular momen- 
tum from it. For equal-mass perturbers, the net drag force in 
the azimuthal direction is thus smaller than the isolated coun- 
terpart. The contribution of the companion wake to the DF 
force in the azimuthal direction is delineated in the inset of 
Figure |2l In the supersonic range, the ratio -X^^^jTi^^ is on 
average ~ 40%, and drops to ~ 12% dX Jv[ ~ 1.2-1 A where 
the net azimuthal drag force is maximized. Since Xi^^ in- 
creases with decreasing ln(rj^ij^/Rp) for A^ > 1, the effect of 
the companion wake on the orbital decay of supersonic dou- 
ble perturbers would diminish as the perturber size nominally 
represented by rmin decreases relative to the orbital radius. In- 
terestingly, in the subsonic case, the ratio -22,(p/Xi,(p, which 
does not depend on rmin, is ~ 50% at M ^ 0.7 and steeply 
approaches unity as A^ ^ 0, suggesting that the effect of the 
companion wake is larger as the speed of perturbers decreases. 

On the other hand, Ii^r and X2^r are of comparable ampli- 
tude over a wide range of M, and thus give rise to a net radial 
drag on double perturbers that is about twice larger than in the 
corresponding single-perturber cases^. They affect the orbital 
eccentricity rather than removing angular momentum much 
(e.g.. Paper I), which may have a gravitational wave signature 
detectable with Laser Interferometer Space Antenna (LISA) 
(Armitage & Natarajan 2005). 

For practical purposes, we fit our results for X2^r and X2^^ 
using 






' 0.5- 0.43 ('l-cosh"^-^^(2.2A^)') if A^<2.97, 
0.76 - 0.08 ^^ + (^-2.76)-^') if 2.97<M<6.2, 

0.56 - 0.027 (m + {M- 6)-!] if M> 6.2, 

(5) 
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Fig. 2. — Gravitational drag forces on a perturber in a double system in the 
radial (blue) and azimuthal (red) directions as functions of the Mach number 
A4. The dashed curves adopted from Paper I with r^m/Rp = 0. 1 represent the 
forces (Z\) originated from the wake of the perturber itself, while those (X2) 
from the companion wake are plotted as dotted lines. The solid lines give the 
net DF forces (Z\ +X2) for equal-mass perturbers. The inset plots the ratio 
-T2^c^/Xi (^ which is positive and less than unity for all A4. 
and 



-0.022 {10- M) tanh(3A^/2) 
-0.13 + 0.07tan-H5A^-15) 



ifA^<2.97, 
ifA^>2.97, 



(6) 



^ When A4 = 0, the steady-state solution of equation ([T) is simply a = 
(GMp/c])(\r-Rp\-^ +f\r+Rp\-^), for which equation (|4) yields Zi^r/M^ = 
0.5mdX2^^/M^=0. 



which are accurate within 6% of the numerical results for all 
M. The algebraic fits to Xi^r and Xi^^ are given by equations 
(13) and (14) of Paper I . It can be shown that Xi^^ ^ ~^2,^ -^ 
A^^/3 in the limit of small M. 

4. DISCUSSION 

For a system composed of two perturbers moving on copla- 
nar circular orbits at the opposite sides of the system center, 
we have found that a perturber is dragged backward by its 
own induced wake, while it is simultaneously pulled forward 
by the wake of its companion. For equal-mass perturbers, the 
ratio of the positive torque from the companion wake to the 
negative torque from its own wake varies between 0.1 to 1.0, 
and has a mean value at about 0.4. This indicates that, since 
the wake tails are a large-scale perturbation, the effect of a 
companion wake on the orbital decay of a double system is by 
no means negligible except perhaps at M ^ 1.2-1 A. When 
the perturbers do not have the same orbital radius, the posi- 
tive torque generated by the companion wake is expected to 
be enhanced (thus making the net azimuthal drag reduced) if 
the companion has a larger orbital radius. In order to quantify 
this effect, we have computed the drag force, for instance, on 
a body with M = 1 ai Rp when the companion is at 2Rp and 
A^ = 2 (in order to have the same Q), and found that the pos- 
itive torque increases by a factor of ~ 1.4. If the perturbers 
have different orbital frequencies so that one perturber com- 
pletes several orbits in the orbital period of the companion, the 
orbit- averaged torque exerted by the wake of the companion is 
likely to be reduced. In the limit of very different frequencies, 
the orbit-averaged torque by the companion wake becomes 
negligible. 

The results of this Letter can b e immediatel y applied to the 
numerical models considered in lEscalaetaP (|T004,) for DF- 
induced mergers of supermassive black holes in a gaseous 
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0.0 0.5 1.0 1.5 2.0 

time [Myr] 
Fig. 3. — Decay of the separation of a black hole binary of equal mass 
caused by dynamical friction due to a background gas. The thick and thin 
curves correspond to the cases when the wakes of both binary components 
are considered and when the companion wake is ignored, respectively. In 
both cases, the solid lines plot the results based on both radial and azimuthal 
forces, while those with only azimuthal forces are given as dashed lines. 

medium. In their model, two black holes of equal mass are 
initially separated from each other widely and orbit at near- 
transonic speed under an external gravitational potential and 
begin to undergo gaseous drag. Our results shown in Fig- 
ure |2] suggest that the companion wake presumably plays a 
minor role in this early phase of the orbital decay. Escala et al. 
(|2QQ4)) reported that at some point when the binary separation 
is reduced to ~ 7 pc, the binary produces a wake in the sur- 
rounding medium that is well approximated by an ellipsoid 
with an axis ratio of 2:1. The major axis of the ellipsoid lags 
behind the binary axis by ~ 22.5°. Since the wake configu- 
ration of this sort can be obtained by blurring the perturbed 
density shown in Figure [T^, we infer that the black holes have 
M ~ 0.6 at this time. The binary keeps hardening by an el- 
lipsoidal torque, and the effect of the companion wake is now 

non-negligible at all. 

Figure 10 of lEscala et al.l (|2004|) shows that it takes the bi- 
nary about 1.5 Myr to decay from 7 to 0.7 pc. To check 



whether this is consistent with our predictions, we consider a 
binary black hole with mass Mp = 5x lO^M© each, embedded 
in a uniform medium with number density ^o = 1 -5 x 10^ cm~^ 
and sound speed Cs = 650 km s~^ , a condition similar to those 
in Esc ala etal.1 (l2004 . Initially, the binary has A^ = 0.6 and 
a separation of ^-q = 7 pc. We follow the orbital decay of the 
binary subject to the DF forces found in ^3.21 Figure [3] plots 
the resulting temporal changes of the separation s. The thick 
curves correspond to the cases when the wakes of both binary 
components are considered, which show that the binary de- 
cays to s = . 7 pc i n ~ 1.5 Myr, consistent with the results of 
lEscalaetaP (|2004)) . When the effects of the companion wake 
are neglected, the decay time becomes shorter by about a fac- 
tor of 1.5, as indicated by the thin lines. Note that the close 
agreement between the two cases with and without the radial 
force demonstrates that it does not affect the decay much. 

The fact that the companion wake always acts against the 
orbital decay of a double system suggests an interesting pos- 
sibility that a less massive perturber in an unequal-mass sys- 
tem may be able to experience a forward net force and move 
radially outward by acquiring (instead of losing) angular mo- 
mentum from the wake of a more massive body, resembling 
a dynamical barrier, while the latter is little influenced by 
the former. This dynamical boost happens if the mass ratio 
/ of two components exceeds l^i (p/22,(p|. The evolution of 
unequal-mass perturbers is more complex since the assump- 
tion of identical orbits for both perturbers will fail soon. In 
practice, one expects that the orbit of the light perturber will 
be quenched to the massive perturber; the former is trapped 
on a radius where the drag of its own wake is slightly larger 
than the forward force of the companion wake, until the orbit 
of the massive perturber decays and the wakes decouple. We 
will discuss this case somewhere else. 
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